Observation of Co-tunneling in Pairs of Coupled Flux Qubits 
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We report measurements of macroscopic resonant tunneling between the two lowest energy states 
of a pair of magnetically coupled rf-SQUID flux qubits. This technique provides a direct means 
of observing two-qubit dynamics and a probe of the environment coupled to the pair of qubits. 
Measurements of the tunneling rate as a function of qubit flux bias show a Gaussian line shape that 
is well matched to theoretical predictions. Moreover, the peak widths indicate that each qubit is 
coupled to a local environment whose fluctuations are uncorrelated with that of the other qubit. 



Superconducting circuits have played an essential role 
in realizing quantum mechanical phenomena in macro- 
scopic systems. One such example is the observation of 
macroscopic resonant tunneling (MRT) of magnetic flux 
between the lowest energy states of single rf-SQUID flux 
qubits, as demonstrated by several groups These 
measurements provide both a clear signature of quantum 
mechanics in a macroscopic circuit at a finite tempera- 
ture and in the presence of noise and a direct means of 
determining the tunneling energy between states. Theo- 
retical descriptions of the MRT rate have been presented 

@ and indicate a direct connection between the pro- 
file of the MRT rate peaks and properties of the envi- 
ronment. Analogous measurements of the tunneling of 
magnetization in nanomagnets 0, 0] suggest that MRT 
is responsible for dynamics in these materials as well. 

In this work, we extend measurements of MRT to in- 
ductively coupled pairs of flux qubits. We present exper- 
imental observations of tunneling between the two lowest 
energy states of the coupled system for several coupling 
strengths. These data yield two-qubit energy gaps that 
match those predicted by the independently calibrated 
Hamiltonian of the coupled system. Moreover, measure- 
ments of the two-qubit energy gap are used to infer sin- 
gle qubit energy gaps at ~ ft. x 10^ Hz without the use 
of microwave lines. Finally, the profile of the MRT rate 
versus qubit flux bias has a Gaussian lineshape with a 
width that is a factor of a/2 larger than that of a single 
qubit. We argue that this observation indicates that the 
environments coupled to each qubit are uncorrelated. 

For a single flux qubit, an MRT experiment consists 
of measuring the rate of tunneling of flux between two 
wells of the double-well potential of the rf SQUID when 
the lowest energy levels of each well are closely aligned. 
Restricting the dynamics of the single rf SQUID to its 
two lowest energy states allows one to map the physics 
of this device onto the canonical qubit Hamiltonian: 
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where ax,z are Pauli matrices, e = 2/p(<i>^ 
energy difference between the two wells, Ip is the ampli- 



$g) is the 



tude of the persistent current in the rf-SQUID loop, A 
is the tunneling energy, and Q is an operator that acts 
on an environment that generates flux noise in the qubit. 
Here, $^ represents the external flux bias applied to the 
rf-SQUID loop and $° is the degeneracy point. Hamil- 
tonian II]) is valid when |e|, A ^ hup, where fiujp is the 
energy spacing to the next excited state of the rf SQUID. 
For a non-Markovian environment , the initial tunnel- 
ing rate from |0) to |1) (eigenstates of a^) versus e has a 
Gaussian profile, as given by Eq. (2) in Ref. 0. 

A natural extension to the single qubit MRT experi- 
ment is to add a second qubit that is inductively coupled 
to a first qubit via a mutual inductance Afoff . The system 
then has the following low energy Hamiltonian: 
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where J = Mcslpilp2 is the coupling energy and all 
qubit-specific quantities are labeled with i e [1,2]. The 
qubits are ferromagnetically (FM) coupled when J < 0. 
For |J| |ei|, Ai, the two lowest energy eigenstates are 
approximately superpositions of the FM ordered states 
1 00) and |11) (eigenstates of ). One can therefore 

write a two-state Hamiltonian to describe the low energy 
dynamics in this subspace: 

H2q « -i[(ei +e2)T. +5tJ + ^(Qi +Q2)t. , (3) 

where t^^z are Pauli matrices in the above 2-dimensional 
subspace and g is the two-qubit energy gap: 



g = xjJ' + i(Ai + A^r + J(Ai - . (4) 

For the regime Ai, A2 <C 2|J|, Eq. dH) simplifies to 
g w A1A2/2J. If 2| J| ^ TiWp, the nearest excited states 
outside of this subspace are formed from the antiferro- 
magnetically ordered states |01) and |10). If 2\J\ > Tiojp, 
additional levels from the two rf SQUIDs need to be in- 
cluded and g must be evaluated numerically. For all mea- 
surements reported herein, 2| J| <^ hujp, thus justifying 
our use of Eqs. Q and (g]). 



By adapting the derivation in Ref. [Bj to the subspace 
spanned by |00) and |11), we derive a functional form for 
the two-qubit co-tunnehng rate from |00) to |11): 
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■ exp 
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where e' = ei + 62, ep,2g and W2q represent the displace- 
ment and width of a Gaussian peak, respectively, and 
rii_>oo(e') = roo-).ii(— e'). We define the noise spectra) 
density S2q{(-^) for the coupled system as follows: 



52,H = J dte^^'{[Qi{t)+Q2it)][Qi{0)+Q2m) , 

(6) 

where (. . .) denotes averaging over all environmenta' 
modes. We use Eq. ([6]) to calculate W2q as in Refs. 0,01: 
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rfie'"*(Q2(t)Qi(0)) (7) 



where W is the width of a single-qubit MRT peak 0, Q . 
If Qi{t) and Q2{t) are uncorrelated the final integrals 
of Eq. ^ win be zero, thus yielding W2q = V^W. As 
in the single qubit case, the fluctuation-dissipation the- 
orem connects the peak width and displacement to the 
temperature of the environment T: = 2Tep^2q 0- 

We performed measurements on a chip that includes 
eight compound-compound Josephson junction (CCJJ) 
rf-SQUID flux qubits llC| with sixteen pair wise tun- 
able coupling elements [ll| . Figure HJa) shows a sim- 
plified schematic of two qubits connected by a coupler. 
For further details on this circuit see Refs. |1CH13| . The 
chip was manufactured on an oxidized Si wafer with 
Nb/Al/Al203/Nb trilayer junctions and four Nb wiring 
layers insulated from one another with planarized high 
density plasma enhanced chemical vapor deposited Si02. 
We mounted the chip in an Al box on the mixing cham- 
ber of a dilution refrigerator. All measurements reported 
herein were performed at T = 21 mK. 

For every coupled pair of qubits in this circuit we had 
independent time- varying control over the CCJJ flux bias 
^ccii ^^'^ the qubit body flux bias [l^. These signals 
were provided by room temperature current sources with 
cold filtering that limited the bandwidth to 5 MHz. The 
couplers provided mutual inductances Mcs between pairs 
of qubit loops and which could be tuned from 1.8 pH to 
—3.0 pH via a static flux bias [llj. We focus the rest 
of the paper on results obtained from a particular pair 
of qubits we call gl and q2, which we isolated from the 
rest of the circuit by setting all but one of the interqubit 
couplers to Afcff = 0. Wc have reproduced these results 
with the other 15 pairs of qubits on this chip. 
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FIG. 1: (a) Schematic showing two compound-compound 
Josephson junction flux qubits and a tunable coupling ele- 
ment. For the experiments reported on herein, we had inde- 
pendent time varying control over "I'ccjj a^nd <l>q for each qubit. 
The coupling strength was tuned with a static flux bias 
(6) Example eigenspectrum for a strongly FM coupled pair of 
flux qubits. The four lowest lying diabats are |00), |01), 
and 1 10) . Anticrossings that give rise to specific resonant tun- 
neling processes are highlighted with dashed ellipses and have 
been denoted as |a/3) — !> |<57). 



We began our experimental investigation of this chip 
by calibrating all on-chip mutual inductances and qubit 
parameters. Reference [l3| describes these calibration 
techniques for an identical chip. We obtained a qubit 
critical current Ic = 3.38 ± 0.01 ^A, a qubit inductance 

26 ± 1 



CCJJ 



Lq = 338 ± 1 pH, a CJJ loop inductance L, 
pH and a qubit capacitance C = 185 ± 5 fF for q\ and 
q2. 

After the parameter calibration noted above, we mea- 
sured Ip and A as a function of *&ccjj > summarized 
in Fig. [21 We measured Jp by using a second qubit as 
a sensor of coupled fiux, as described in Ref. [10|. We 
measured A via three methods. The first method used 
MRT between the ground and first excited states of a sin- 
gle qubit . Figure [3] shows example single qubit MRT 
rate measurements for three values of ■ We obtained 
W/kB = 26 ± 2 mK {W/2Ip = 80 ± 6 n<^o) from fitting 
these data to Eq. (2) of [2]. The second method involved 
Landau-Zcner (LZ) rate measurements [l^ . The band- 
width restrictions of our cold filtering limited both of the 
above techniques to measuring A//i < 50 MHz. To char- 
acterize larger A, we used an alternate method that will 
be described below. The solid curves in Fig.[2]are the the- 
oretical predictions of the ideal CCJJ rf-SQUID Hamilto- 
nian given in Ref. [loj using the independently calibrated 
qubit parameters Ic, Lq, iccjj and C cited above. The 
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FIG. 2: (a) Ip versus 'i&ccjj for and q2. (b) A versus <I>|f(,jj 
for ql and q2. A was obtained from single-qubit LZ (IQ-LZ), 
single-qubit MRT (IQ-MRT), and two-qubit MRT (2Q-MRT) 
measurements. 



quality of the agreement between experimental data and 
theory justifies our identification of these devices as flux 
qubits. 

With the single qubit calibrations completed, we then 
turned to experiments on pairs of strongly FM coupled 
qubits. The two-qubit MRT experiment was performed 
in a manner very similar to the single qubit case 0. We 
applied ±10 m$o to each qubit in the presence 

of ^ccj) = — *^o/2, raised their tunnel barriers by ramp- 
ing from — $o/2 to — and then waited 1 ms to 
ensure the two-qubit pair was in its ground state. This 
initialized the coupled pair in either |00) or with cer- 
tainty. Next we adjusted ^ and then simultaneously 
lowered the tunnel barriers of both qubits for a dwell 
time T before again raising them via the individual $ccjj 
. We measured the loss of probability from the initial 
state and repeated for a range of t. The probability of 
the initial state as a function of r was fit to an expo- 
nential to extract roo^-ii or Fn-^oo depending upon the 
initialization. 

Having individual control of <&ccjj ^^^r each member of a 
pair allowed us to perform two-qubit MRT measurements 
in which we either matched Ai = A2 or deliberately 
mismatched Ai 7^ A2. For the first set of measurements, 
we set Ai = A2 by biasing ^'ccjj.i = '^'?cjj.2- Figure 
m shows example measurements of Foo^ii and Fn^oo 
as a function of $^ ^ with = and Mcs = -2.35 
pH. The increase in rate for l^^^l > 0-4 m$o is due to 
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FIG. 3: Example measurements of single qubit MRT rates 
versus for q2. The hollow (solid) symbols are Fo^.! 
(Fi^o). Data shown are for $?cjj/'5o = -0.6231,-0.6242, 
and —0.6253 from top to bottom, respectively. The curves 
are fits to Eq. (2) of Q]. 



H^2g/fcs = 39±2mK 
r = 22 ± 2 inK 
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FIG. 4: Example measurements of two-qubit MRT rates ver- 
sus ^q^i) for two coupled qubits (Mos = —2.35 pH). The 
hollow (solid) symbols are Foo^ii (Fn^oo). Data are shown 
for both qubits biased at 3>ccjj/<E>o = -0.6176, -0.6181, and 
—0.6187 from top to bottom, respectively. The curves are fits 
to Eq. ©. 

tunneling from the initial state, either |00) or |11), to 
|01) or |10), by the processes depicted in Fig. [Tfb). For 
these experimental settings, it was predicted that these 
processes would peak at ±0.8 m<i>o. 

To extract ep,2q and W2q from data such as those in 
Fig. m we fit the MRT rate peaks to Eq. (P. For all 
measurements with Ai = A2 we obtained 14^213/^5 = 
39 ± 2 mK. The ratio W2q/W = 1.5 w which 
indicates that the environment coupled to gl is uncorre- 
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FIG. 5: Two qubit gap g versus A for a range of Mcs- The 
qubits were biased at matching 'li'ccjj • Values of g were ob- 
tained from fits such as those shown in Fig. |4] The sohd 
curves are the theoretical predictions of from Eq. Q. 

lated with that coupled to q2. This is evidence that the 
source of flux noise in these qubits is local to the qubit 
wiring, which agrees with the conclusions of others [l5l |. 
Values of W2q and ep^2q were used to infer a temperature 
T = W^^/2ep^2q = 22 ± 2 mK. This is consistent 
with T as determined via single qubit MKT and with 
that reported by thermometry. 

Fitting MKT peaks to Eq. [5] also allowed us to extract 
the two-qubit energy gap g. Figure [5] shows g for a range 
of single qubit A = Ai = A2 and three different cou- 
pling strengths. The theoretical predictions were gener- 
ated using Eq. ([4]), where we used A as predicted by the 
theoretical curve shown in Fig. [21(b) and an independent 
calibration of J. There is good agreement between the 
measured and predicted g for different Moff. We conclude 
that (7 (X A^ and g oc 1/J, as predicted when A ^ 2| J|. 

Besides having achieved the goal of demonstrating 
quantum mechanical behavior in pairs of coupled qubits, 
we have used two-qubit MKT measurements to comple- 
ment our single qubit methods for calibrating large A. 
We biased control qubit ql at $|?j,jj/$o = -0.6222, a 
point at which Ai/h = 6.0 ± 0.2 MHz had been in- 
dependently calibrated using single qubit methods. We 
then targeted ^^^^^ /$o > -0.6222 of the second qubit 
q2, which ensured that Ai < A2. The dynamics of the 
coupled system were then governed hy g < Ai, which 
made them measurcablc given our experimental band- 
width. We extracted g from a two-qubit MKT measure- 
ment of this mismatched pair which, along with the in- 
dependently calibrated values of Ai and J. allowed us to 
infer A2 from Eq. Q. We then switched the roles of ql 



and q2 to infer Ai. The results of such experiments us- 
ing Moff = —2.91 pH have been summarized in Fig.[2l^b). 
Again, the experimental data agree with the predictions 
of the independently calibrated ideal CCJJ rf-SQUID 
Hamiltonian, further confirming the self-consistency of 
our measurements. Thus we have validated a new tech- 
nique for characterizing high tunneling energies of single 
qubits despite the limited bandwidth of our apparatus. 

Conclusions. Macroscopic resonant tunneling is a pow- 
erful way of characterizing single and coupled pairs of 
superconducting flux qubits. We have demonstrated 
that inductively coupled pairs of flux qubits behave as 
expected by quantum mechanics in that the two-qubit 
tunneling energy g inferred from fitting experimental 
data agrees with the predictions of an effective two-level 
Hamiltonian for the coupled qubit system. Investigat- 
ing one pair in detail for a range of Afos, Ai, and A2 
yielded g that not only matched theoretical predictions, 
but allowed us to probe single qubit tunneling energies 
0.2 GHz < A/ft, < 2 GHz without the use of microwave 
lines. Finally, the widths of the two-qubit MKT rate 
peaks were a factor of larger than that of a single 
qubit. It was argued that this is an indication that the 
environment interacting with one qubit is uncorrelatcd 
with that of the other qubit. This latter observation im- 
plies that the source of flux noise is local to the qubit. 
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